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Abstract 



o 

C3 . tional correlation by random but large amplitude molecular jumps, as opposed to infinitesimal small 



Ivanov- Anderson (IA) model (and an earlier treatment by Kubo) envisages a decay of the orienta- 



i ' jumps assumed in Brownian diffusion. Recent computer simulation studies on water and supercooled 
£h ' liquids have shown that large amplitude motions may indeed be more of a rule than exception. Ex- 

o : 

, O, - isting theoretical studies on jump diffusion mostly assume an exponential (Poissonian) waiting time 
, distribution for jumps, thereby again leading to an exponential decay. Here we extend the existing for- 
f^. ' malism of Ivanov and Anderson to include an algebraic waiting time distribution between two jumps. 
t-h ■ As a result, the first {I = 1) and second (I = 2) rank orientational time correlation functions show 

cn ; 

^5 . the same long time power law, but their short time decay behavior is quite different. The predicted 
I — ' Cole-Cole plot of dielectric relaxation reproduces various features of non-Debye behaviour observed 

o : 

' experimentally. We also developed a theory where both unrestricted small jumps and large angular 
^ ■ jumps coexist simultaneously. The small jumps are shown to have a large effect on the long time 

b : 

5t , decay, particularly in mitigating the effects of algebraic waiting time distribution, and in giving rise to 
an exponential-like decay, with a time constant, surprisingly, less than the time constant that arises 
from small amplitude decay alone. 

PACS numbers: 05.40+jc, 64.70.Pf, 66.20.+d, 66.10.-x 
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1. INTRODUCTION 



Because orientational relaxation of molecules is relatively easily accessed by a variety of 
experimental techniques (NMR, IR, fluorescence depolarization, optical Kerr effect, to name a 
few) [3, 3), 0, [5, Q], many theoretical models and microscopic studies have addressed various 
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id, 



111 ] . The most celebrated 



aspects of molecular rotation in dense liquids and glasses 
of these studies is the one carried out by Debye many years ago, in terms of a simple rotational 
diffusion equation. 0, 4] The theory assumes that rotational correlation in liquids decay by small 
amplitude rotational Brownian motion, with a rotational diffusion coefficient, .Dr. The theory 
makes simple prediction that the decay of the correlation functions of all ranks of Legendre 
function is exponential and is given by 



C t (t) = {P £ [cos0(t)]> = exp{-£{£ + l)D R t). 



The Debye model of rotational diffusion has played a pivotal role in most of the discussions 
on rotational diffusion in the condensed phases. Many theoretical studies have attempted to 
improve upon Debye model, using for example, a generalized Langevin equation approach which 
ultimately leads to a time or frequency dependent diffusion coefficient. But all these approaches 
use the general assumption of infinitesimal rotation of Brownian motion. 

Recent experimental and theoretical studies on supercooled liquids and surpri sing ly, liquid 



water, have shown a marked departure from the classical Debye behaviour. 
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131 ] In these 



cases, rotational diffusion is found to occur by large amplitude jumps. In supercooled liquids 
and glasses, the orientational relaxation is often markedly non-exponential. 12] In liquid water, 
the non-exponentiality, at room temperature, is weak but relaxation becomes progressively 
non-exponential at low temperatures. [3] The rotational relaxation (and also translational 
diffusion) seems to occur primarily through rare but large amplitude jumps. A quantitative 
understanding of the origin of such jumps has remained a subject of great interest, though 
largely unsolved. 

A model of rotational jump diffusion was actually proposed by several workers in the past, 
most notably by Kubo 14, 15] and Ivanov Q. In Kubo's model of jump diffusion, the rotator 
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was restricted to jump in a circle, that is restricted to two dimension. (See Figure 1(a)) Ivanov's 
model was more general where jumps were isotropically distributed in three dimension. (See 
Figure 1(b)) Ivanov's model has found increasing use in describing experimental results. In 
this model, the waiting time between jumps obey an exponential distribution and as a result, 
the decay is single exponential, with the first and second rank correlation functions which are 
given by 



where A is the constant amplitude of jump and r is the average time interval between any two 
jumps. The probability that there are n jumps in a time interval t is assumed to be given by 
Poisson distribution 



For small jumps, equations 2 and 3 go over to the Debye behaviour, with ri/r 2 = 3, where T\ and 
r 2 are the decay time constants of Ci(t) and C 2 (t), respectively. However, the relaxation pattern 
is different for long jumps. The difference is most acute when the jump angle A is close to 7r. A 
jump by A ~ 7r relaxes C\{t) but not C 2 (t). therefore, T1/T2 becomes much smaller than 3 and 
approaches zero. Interestingly, for intermediate values of jump length parameter (A « 7r/2) , the 
ratio of the two time constants approach unity. Simulations have often found that the values of 
the ratio ri/r 2 lies close to unity in supercooled liquids which has been taken as an indication of 
the emergence of jump diffusion as a contributing mechanism of rotational relaxation. On the 
other hand, in a notable development, recently Dejardin and Jadzyn have extended the Debye 
rotational diffusion model to fractional rotational diffusion case to treat events that are non- 
local in time due to the existence of extensive memory effects. [16| Under fractional rotational 
Brownian motion, the decay of the relevant orientational correlation functions is no longer 
exponential (or, sum of a few exponential terms). It was shown that the decay can be described 
as a combination of Mittag-Leffler temporal pattern, behaving like a stretched exponential at 
short times and an inverse power law in the long time limit. Fractional equations provide an 
anomalously slow decay at long times, often referred to as subdiffusive regime. As pointed 
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(3) 




(4) 
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out by many, fractional Fokker-Planck or Smoluchowski equation is a natura 



generalization of 



normal diffusion to disordered systems with scale free memory effects. |17|, ll8|, ll9|, |20] Dejardin 
and Jadzyn obtained analytic expressions of the frequency (a?) dependent electric birefringence 
spectrum, x{ u )- They plotted Cole-Cole diagram for various cases to demonstrate the effects 
of fractional diffusion. Fractional diffusion gives rise to markedly non-Debye Cole-Cole plot 
which now varies from Cole-Davidson skewed arc behaviour to Cole-Cole depressed circle. The 
treatment of Dejardin and Jadzyn is still based on rotational diffusion equation, and does not 
contain the effects of finite jumps. However, as already discussed, in many complex systems, 
rotational relaxation occurs by large amplitude jumps. For example, in liquid water, it has 
recently been discussed how much of the relaxation of orientational correlation occurs by jumps 



which are of the order of 60°. 



13] Such motions cannot be treated as a Brownian diffusion. In 



addition, in many complex systems, the waiting time distribution between jumps may not be 



approximated by Poissonian distribution. [12|, 
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22] In the present work we have 



extended the theory of Kubo, Ivanov and Anderson to treat jump diffusion with an algebraic 
waiting time distribution. Since the treatments of Kubo are different from that of Ivanov and 
Anderson, separate solutions have been obtained. As is usually the case for relaxation with 
fractional diffusion, an analytical solution of the orientational time correlation function has 
been obtained only in the frequency domain. However, we have been able to obtain asymptotic 
solution in all the cases. The decay of the correlation of the first and the second rank harmonic 
follow power law. This signifies a breakdown of Debye behaviour in dispersive medium. We 
have also developed a theory where the large amplitude jumps simultaneously coexist with 
small amplitude jumps. Interestingly, the exponential kinetics is recovered even under the 
small fraction of small amplitude jumps but the decay is accelerated by large amplitude jumps. 
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2. GENERALIZATION OF JUMP MODEL 



A. Kubo model 

We consider a 2-dimensional rotator which makes the series of jumps. The model was 
originally proposed by Kubo for rotational relaxation of spin under the presence of pulsed 
magnetic fields, [ijl [ij]] The angle changes by the jump with the amount, A n . According to 
Kubo we define, 

f(t) = /f>xp(i f> m )P(M)\, (5) 

\n=0 \ m=0 / / 

where P(n, t) represents the probability of having n jumps up to time t and (• • • ) denotes the 
ensemble average of scattering angles. The correlation functions are obtained by, 

Ct(t) = Re [fit)] , (6) 

where £=1,2. It should be noticed here that the I— th rank correlation function is defined in 
terms of cos £8(t) for two dimensional rotator. Even when the second rank correlation function 
is defined in terms of Legendre function of three dimensional isotropic rotator, the right hand 
side of eq. (ED is equal to (Ci(t) — Cg(oo)) / (Cg(0) — Cg(oo)), which still expresses the decay 
of correlation function. We calculate Re and represents it by Cz(t) of 2 dimensional 

rotator for Kubo model. The waiting time distribution of each time interval r n = t n — £ n _i, 
is assumed to be statistically independent and it is represented by, ipij). We also introduce 
(p(t) = J t °° dtiip(ti) as a probability that a molecular rotor will not make a jump for the time 
interval between and t. The probability of having n jumps up to time t is given after the 



Laplace transformation, P(n,s) = J °° exp(—st)P(n,t), as [23] 

P(n,s)=0(s)r(s), (7) 

where <p(s) = (1 — ip(s))/s. The amount of each jump is also assumed to be statistically 
independent, 

(exp(z \J A n )) = exp(mA). (8) 



Then, we have, 



oo 



C e (t) = Re ^2exp(i£nA)P(n,t) 



(9) 



n=0 



By substituting eq. ((Tj) into eq. dHJ), the Laplace transform of it becomes, 




l-V(s) 



Re 



1 



(10) 



1 - exp(z£A)^(s) 



The above equation holds for an arbitrary waiting time distribution function and an arbitrary 
jump amplitude. When the waiting time distribution decays exponentially with the charac- 
teristic time of the jump interval, r, the well known result for the Poisson noise is recovered, 
i.e., 



Eq. (TlOT) is the generalization of Kubo model to the case of an arbitrary waiting time distri- 
bution. Eq. pip was given by Kubo for the Poissonian waiting time distribution for £ — 1. 
"bl flo1 | In 2 dimensional rotational diffusion, eq. (p]) is replaced by exp (—£ 2 D R t) j^J since the 
correlation function is defined by (cos (£#(£))). Eq. ( TTTT) gives T1/T2 = 4 in the limit of small 
amplitude jump, which is consistent with the result of conventional 2-dimensional rotator. 

B. Ivanov- Anderson model 

For dielectric relaxation, however, Ivanov-Anderson model is more popular to describe the 
random scattering of angle. Ivanov-Anderson model is, however, similar in spirit to Kubo 

model. Ivanov-Anderson model assumes for simplicity an isotropic reorientation by random 
angular jumps. If we denote the angular change by n-th jump by A n and assume that it is 
statistically independent as in Kubo model, we have, 



C e (t) = exp [-(t/r) (1 - cos(M))] cos [(t/r) sm(£A)} . 



(11) 



L\ = cos A = (cosA n ) ; 

L 2 = (3 cos 2 A - l)/2 = ([(3 cos 2 A n - l)] /2). 



(12) 



(13) 
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It is known that the correlation function is expressed in the form similar to eq. as, 

oo 

C e (t) = J2^P(n,t). (14) 

n=0 

By substituting the Laplace transform of eq. ([7]), the Laplace transform of the correlation 
function, eq. ( fT4l) . becomes, 

c e (s) = — (15) 

si- Ld>{8) 

Eq. (1151) is valid for an arbitrary waiting time distribution. When the waiting time distribution 
between jumps decays exponentially with the characteristic time r, eq. f[T5"j) recovers the well 
known results of eqs. (J2J)- (EI) - Here, the Ivanov- Anderson model is generalized to the case of an 
arbitrary waiting time distribution. 



3. ALGEBRAIC WAITING TIME DISTRIBUTION OF JUMP 

Now, we study the influence of the power law waiting time distribution ip(t) ~ l/t a+1 on 
the correlation functions. A well-known and popular waiting time distribution function which 
is normalizable is given by, 

m - a7( ;;i; 7 ' t) . (») 

/r 



where 7(z,p) = Jq e~ l t z ~ l dt for (Re^ > 0) is the incomplete Gamma function. [25j Eq. (|T6l) 
can be derived from a model in which the hopping rate depends on a parameter exponentially 
for example activation energy) and the value of this parameter has an exponential distribution. 



261 ] 7 r is a parameter characterizing the hopping frequency. The waiting time distribution of 



eq. (I16p can be written approximately in a more transparent form, 



, , . aT (a + 1) 

7r 



where T(z) is the Gamma function, 
is obtained as, 



25| The Laplace transform of the waiting time distribution 



i>(s) = 1 - 2 Fi [l,a,a + 1, -j T /s] 



Eq. (jTSJ) is useful because the following property of the hypergeometric function is known [25] , 



^[1,0,0 + 1,-7,/*] «-^( -J -<1- (19) 

sm7ra \7r/ 7r 

By substituting eq. ffTSj) into eq. ffTUl) . extended Kubo model leads to, 

C t (s) ~ 2F1 [i,a,q + 1,-7^]^ / 1 \ (2Q) 



7TQ; 

sin 7TQ; 



1 - cos(£A) 



1 - exp(i£A) 
1 



s l-cw* ■ 



(21) 



^l-cos(IA)) z + sin i (^A). 
By performing the inverse Laplace transformation, we find the asymptotic time dependence for 
extended Kubo model, 



sin na 



1 - cos(£A) 



' (22) 



1 - cos(£A)) 2 + sin 2 (£A)J T(l - a) { lv t) a ' 
We can obtain the solution of the extended Ivanov- Anderson model in a similar fashion. By 
substituting eq. (fl9l) into eq. (fTB"l) . eq. (fl5]) becomes, 



In the long time limit, 



1 s Q_1 

<^(s) = t = J—. . (23) 

L/£ n 1 Lg sin 7rcx K ' 

S + —j 7r 



1 — Le sin to 



we find, in the small s limit, the following asymptotic result, 

Ct{s) ~ * -J_. (25) 

By performing the inverse Laplace transformation, we find the asymptotic time dependence of 
extended Ivanov- Anderson model as, 

Ct{t) « ™ h — (26) 

sm7ra (1 — L e ) T{1 — a) (jj-t) 

Both Kubo and Ivanov-Anderson models predict exactly the same time dependence of the 
correlation functions, which is described by the algebraic time dependence with the exponent 
a. As we can see from eq. ( )24l) . the above asymptotic dependence is hardly attained when 
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Li w 1. This case can be examined more rigorously. By the inverse Laplace transformation of 
eq. (123j) . we find the correlation functions expressed in terms of Mittag-Leflier function, 

1 — Lt sin ira 



C t (t) = E a 



drty 



(27) 



where the Mittag-Leffler function is defined by, E a (z) = z k /T(ak + 1) and Le ~ 1 is 
introduced. Mittag-Leffler function is approximated by the stretched exponential function, 

1 — Le sin(7ra) 



Ct(t) « exp 



drty 



28 



T(l + a)Li na 

except the final component showing algebraic decay of eq. (j26|) which appears at very long 
times when L< Ri 1. Li ~ 1 occurs when A ~ for both £ = 1 and £ = 2. When A ~ 7r, 
L2 ~ 1 but Li < 1. Therefore, when A is small, both Ci(t) and C^i) mainly decay according 
to stretched exponential law. As A is increased into the range, A 7r, both C\(t) and 
C^if) exhibit algebraic decay of eq. (|26|) . When A is closed to ir, C%(t) mainly decays by the 
stretched exponential law and the final small components of both C\(t) and C*2(t) are described 
by algebraic decay. 

By the known identity of the Mittag-Leflier function for a = 0.5, eq. ( 1271) can be expressed 

as, 

2 



C e (t) = exp 



2(1- L, 



irLf 



erfc 



2(1 



TlLf 



(29) 



for a = 1/2, where the complementary error function is defined by erfc(z) = 
(2/^/W) f°° exp(— y 2 )dy. [25] We compare these analytical results with the numerical Laplace 
inversion of the exact results (eqs. ( TTUT) and ( |T5l) with eq. ( fl8i) ) in sectional 
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4. SIMULTANEOUS COEXISTENCE OF LARGE AND SMALL AMPLITUDE 
JUMPS 



In many situations, jumps with large amplitudes are rare and they are superimposed on 
the small amplitude jumps which occur frequently. We denote the waiting time distribution of 
large amplitude jump by, ip a if) and that of small amplitude jump by ip b {t). Theoretically, they 
are related with each other. By denoting the waiting time distribution of large amplitude jump 
in the absence of small amplitude jump by, ifj^\t) } and that of small amplitude jump in the 
absence of large amplitude jump by, ifij® (t) , the waiting time distribution under the presence 



of both types of jumps is expressed as, 



23| 



Mt) = ^a ] (t) ( 1 



n,(t) ^ 0) (*) ( 1 - 7 ^i 0j (*))- 



dt4°\t) 



(30) 
(31) 



since the large amplitude jump occurs before the occurrence of small amplitude jump and vice 
versa. 

First we consider the case of Kubo model. The Laplace transform of the correlation function 
is generalized from eq. ffTUj) as, 



CAs) 



1 - $a(s) - 4>b(s) 



1 - 4>a{s) - 4>b{s) 



Re 



Re 



^2 ( ex P (^A a ) 4> a (s) + exp (i£A b ) ip b (s 
1 



ra=0 



(32) 
(33) 



1 - exp (i£A a ) ip a (s) - exp (i£A b ) ip b (s) 

where we assume that each jump belonging to either large or small amplitude is statistically 
independent and the average is denoted by, exp (i£Aj), where j = a, b. 

For Ivanov-Anderson model, the Laplace transform of the correlation function is generalized 
from eq. ( 1T5|) as, 

j) a (s) - tp b (s) 
s 

_n=0 

1 



CAs) 



1p a (s) - ljj b (s) 



1-L^ a (s)-L^ b (s) 



(34) 
(35) 
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where we again assume that each jump belonging to either large or small amplitude is statisti- 
cally independent and the average is denoted by, 

L[ j) = cosA i; (36) 
hf = (3cos 2 A.,--l)/2, (37) 

where j = a, b. 

When the both amplitude jumps occur according to the exponential waiting time distribu- 
tion, 



ip a (t) = — exp (-t/r a - t/r b ) 
tp b (t) = — exp (-£/ r a - t/ n) 

n 



Kubo's result of eq. (TTTj) is generalized to, 



C e (t) = exp 



£ (t/rj) (1 - cos^A,-)) 



J=a,i 



COS 



sin^A,) 



J=a,t 



(38) 
(39) 

(40) 



While, Ivanov-Anderson's results, eqs. 



are generalized to, 



C\{t) = exp 



C 2 (t) = exp 



- (l-cosA,-)^ 

j=a,b 



£§(» 



COS 



j=a,b 



(41) 
(42) 



as expected from the Markovian kinetics. 

When, the large amplitude jump obeys the algebraic waiting time distribution, 



cry (a + l,t/r a ) 



t(t/r a y ' 

and the small amplitude jump obeys the exponential kinetics, 

4°\t) = Iexp(-t/r 6 ), 

n 

the waiting time distributions under the presence of both types of jumps are given by, 

Mt) = ^\t)exp(-t/r b ), 

1 f°° 
Mt) = -exp(-t/r 6 ) / dt^Xh). 
n J t 



(43) 

(44) 

(45) 
(46) 
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The Laplace transform is obtained as, 



M S ) 

Ms) 



^\s + l/r b ) 



1/t 6 



i/n 



(47) 
(48) 



where ^^(z) is expressed as, 



^\z) = 1- 2 F 1 [1, a, a + 1,-1/ (r a z) 



(49) 



By substituting eqs. fj4"T|) and (jl8]l into eqs. fl33l) and fl35l) . the exact solutions are obtained in 
the Laplace domain. However, the inverse Laplace transform is too complicated for analytical 
calculation and will be performed numerically. 

Since we are interested in the case where the characteristic time of small amplitude jump is 
much shorter than that of large amplitude jump, the approximate solutions can be obtained by 
taking the limit, 2 F\ [1> a, a + 1, —1/ (t (s + l/r h ))} « 2 F\ [1, a, a + 1, — rj/rj, which leads to 
"4>a{s + 1/ti) ~ ■^i°' ) (l/r fe ). In the approximation, the s-dependence in the Laplace transform 
of waiting time distribution is dropped out. The s-dependence can be ignored when we are 
interested in the time scale larger than that of small amplitude jumps. The quantity is just 
the time integration over waiting time distribution of large amplitude jumps and represents 
the probability of occurrence of large amplitude jumps, i.e. the escape probability from small 
amplitude jumps by a large amplitude jump. By introducing the approximation, we find that 
the correlation functions decay exponentially, 



Ci(t) 

c 2 (t) 



exp 



exp 



(t/n) l 



-3 (t/ n ) 



cos (At) MO) 
1-cos (A a )4(0) 

cos(2A a )^ a (0) - cos(2A 6 )^(0) 



(50) 
(51) 



3 + ^ b (0) -3cos(2A a )^ a (0) 

When A a is close to At ~ 0, we recover eqs. ([2J-© with A = and r = by substituting 
^&(0) w 1 and ^a(O) ~ into the above expressions. The exponential kinetics of eqs. (1501) - ( 15T1) 
even under the algebraic waiting time distribution of the rare but large amplitude jumps is the 
important result of this article. 
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5. NUMERICAL RESULTS 



Fig. [2] shows the relaxation of correlations of Ivanov- Anderson model for A = 60°. For this 
value of A, C^t) decays faster than C\(t) for any values of a. Compared to the decay by normal 
diffusion, decay of both C^t) and Ci(t) is slowed down when the waiting time distribution has 
an algebraic time dependence. At long times, all decay curves are well represented by algebraic 
decay of eq. f)26p with the exponent a. The long time tail in the decay of correlation originates 
from the algebraic time dependence of the waiting time distribution. The exponent is the same 
for both 1 = 1 and I = 1 components. 

When A is increased to A = 170°, C*2(t) decays slower than C\(t) for any values of a as 
shown in Fig. [31 As a is decreased, the relaxation becomes slower. Long time asymptotic decay 
is again represented by the algebraic decay of eq. (|26l) . C 2 {t) decays according to Mittag-Leffler 
function at intermediate times and the fitting becomes better as a is lowered. Since the Mittag- 
Leffler function is approximated by the stretched exponential decay of eq. (1281) in this time 
range, C^if) is well fitted by the stretched exponential function. The stretched exponential 
decay appears when A is close to 180° for C 2 (t) and A is close to 0° for both C 2 (t) and Ci(t). 
For other values of A the analytical form of the decay is not obtained but as we can see from 
Fig. H]the decay curves are very different from the single exponential decay when a < 1. Only 
when a = 1, the correlations decay fast with exponential time dependence which is shown by 
the straight line in the log-linear plot. 

The results for different values of A with a kept constant are summarized in Fig. [51 For 
any value of A asymptotic time dependence is described by the algebraic time dependence with 
the exponent a. The relaxation of Ci(t) slows down monotonically as A is decreased. Since 
C*2(t) is given by the square of cos A, the decay is the same when A is changed to it — A. As 
A is increased from 0°, the decay of C2(t) becomes faster. The fastest decay is obtained for 
A = 90°. When A is further increased from A = 90°, the decay is slowed down. 

In order to facilitate comparison with dielectric relaxation experiments in dispersive medium, 
we investigated the behavior of the Cole-Cole diagram as a function of system parameters (like 
jump angle and value of the power law exponent). The Cole-Cole diagram is obtained from the 
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complex dielectric function, e(o>), which satisfies, [15l. |28| 

e °° = 1 + iwCt (-ico) , (52) 

where e s t is the static dielectric constant and eoo is the optical dielectric constant. The semi- 
circle of Debye relaxation is obtained for normal diffusion. As a is decreased, the maximum of 
e" is decreased and the Cole-Cole diagram is depressed. When A is relatively small, the Cole- 
Cole diagram is symmetric, as shown in Fig. O for any value of a. Such change of Cole-Cole 



diagram by decreasing a has been a 



influence of the jump amplitude. [29|, 1301 ] 



ready noticed, but without paying much attention to the 



Now, we investigate the influence of the jump amplitude on the Cole-Cole diagram. For 
the same value of a, the Cole-Cole diagram is skewed and becomes asymmetric as the jump 
amplitude A is increased as shown in Figs. (7J1HJ The asymmetry is larger for smaller values of 
a. 

The results of Kubo model are quite similar to those of Ivanov-Anderson model. In Kubo 



3raic or stretched exponential 
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15] it has been pointed out 



model, however, an oscillation is observed, in addition to the a, 
decay, as shown in Figs. [9lfT0l In Kubo model of oscillator 
that the jump contributes both to the damping and the oscillation of correlation functions. 
Dipoles rotate in a direction with an average amplitude A, which gives rise to oscillation in 
correlation functions. In the case of magnetic resonance as originally studied by Kubo, the 
oscillation corresponds to the shift of resonant frequency found in the more elaborate theory. 
In the case of electric dipoles, oscillation is due to rather unphysical modeling of rotator, which 
has a preferable direction of rotation. 

The oscillation disappears by decreasing a values. When the jump amplitude is small, Ci(t) 
oscillates more than C^t). For A = 60°, C\(t) oscillates even when a is decreased to a = 0.8. 
In the case of A = 170°, C^t) oscillates for normal diffusion while other decay curves exhibit 
essentially monotonic decay. 

So far, we have investigated the effect of long algebraic waiting time on the rotational 
relaxation with large jump amplitude. Large amplitude jumps are normally superimposed by 
small amplitude jumps. Therefore, we study the case when both large and small amplitude 
jumps are simultaneously present. As shown in Fig. [HJ in this case the kinetics is found to 
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be almost exponential and the decay is well represented by eq. ( 1501) for C\{t) and eq. (15T1) for 
C^it)- The exponential kinetics is also confirmed from the Cole-Cole diagram shown in Fig. 1121 
where semi-circle of Debye relaxation can be seen. The decay curves obtained by multiplying 
eqs. (J2|)- (J3j) for small amplitude jump with the inverse Laplace transform of eq. ([15]) . where 
eq. (1151) for large amplitude jump is introduced, are different from the exact numerical results. 
The exponential kinetics is of course different from the algebraic decay or stretched exponential 
decay found for the large amplitude jump alone. The small amplitude jumps are more frequent 
than the large amplitude jumps. The former gives rise to an exponential decay but in the 
presence of large amplitude jump the correlation functions decay almost exponentially with 
the time constant much smaller than that for the small amplitudes alone. Thus, although 
the large amplitude jumps are rare, they contribute a lot to the relaxation. The waiting time 
distribution of large amplitude jumps in the absence of small amplitude jumps has an algebraic 
asymptotic time dependence which is easily interrupted by the more frequent small amplitude 
jumps. The decay is exponential as a result of cumulative small amplitude jumps, but the 
decay is accelerated by the large amplitude jumps. 

As explained, the exponential kinetics of correlation functions, C\(t) = exp (— t/rx) and 
C%{t) = exp (— t/r 2 ), results from the interplay between the frequent small amplitude jumps and 
the rare events of large amplitude jumps having algebraic waiting time distribution. Although 
the results are not shown, substantially the same results are obtained even when a is lowered 
to 1/2. For a < 1, Tx and r 2 are obtained from eq. (1501) and eq. (15T1) . respectively, while the 
results of normal diffusion are given by, eqs. pTl) - (j42l) . 

We now investigate the jump amplitude dependence of the ratio of the correlation times, 
tx/t2- As shown in Fig. [T3], the result for a = 0.8 and that for normal diffusion are almost the 
same, although the apparent functional forms are seemingly very different. When A a is small, 
the ratio is close to 3 and it decreases monotonically to zero as A a is increased close to 180°. 
The small deviation between the two lines is increased as A a is increased. Since the results 
among different values of a are very close, it could be difficult to judge from the correlation 
functions whether the long amplitude jumps have algebraic long time tail in the waiting time 
distribution. 
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The exponential kinetics is robust as shown in Fig. [TH Even when the frequency of the 
small amplitude jumps is the same as that of the large amplitude jumps, the decay is still 
exponential. However, when the small amplitude jumps occur 10 times less frequently than 
the large amplitude jumps, then non-exponential kinetics is observed in the short time regime, 
t < n . 

6. CONCLUSIONS 

The models of Ivanov and of Kubo are well-known models of the decay of the orientational 
correlation by random but large amplitude molecular jumps. These models are expected to 
give rise to decay dynamics quite distinct from the Brownian diffusion by infinitesimal small 
jumps. When the waiting time distribution for large jumps is Poissonian, even the models of 
Ivanov and Kubo lead to exponential decay, making the detection of appropriate molecular 
mechanism by experimental measurements rather difficult. 

However, if the waiting time distribution is not Poissonian, then the difference between large 
jumps and the Brownian diffusion can be significant. In the present work, we have extended 
the models of Ivanov and Kubo to include such cases. In particular, we have employed an 
algebraic waiting time distribution for large jumps, as such distributions can be useful to 
describe supercooled liquids and glasses, and also in restricted geometries. Recent computer 
simulation studies on water and supercooled liquids have indeed shown that large amplitude 
motions may be more of a rule than exception. 

In the present study, we have solved the theoretical models with algebraic waiting time 
distribution analytically to obtain the first (£ = 1) and second [I = 2) rank orientational time 
correlation functions. As expected, the decay is non-exponential, with power law at longer 
times. The correlation functions for i = 1 and i = 2 show the same long time power law 
exponents, but the short time decay behavior is quite different for the two correlation functions. 

In order to facilitate comparison with experiments on dielectric relaxation, we have calculated 
Cole-Cole plots generated by the orientational decay for a wide variety of parameters, such 
as jump amplitude and the power law exponent. The predicted Cole-Cole plot of dielectric 
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relaxation reproduces various features of non-Debye behaviour observed experimentally. 

In addition, We have developed a theory where both unrestricted small jumps and large 
amplitude jumps coexist simultaneously. The small amplitude jumps are shown to have a 
large effect on the long time decay, particularly in mitigating the effects of algebraic waiting 
time distribution. Thus, in the limit of appreciable number of small jumps, we find the decay 
to become single exponential. We find some what surprisingly that the decay time of this 
exponential decay is much smaller than that given by small jumps alone. 

In this work, we have not made any specific application of our theory but shown that several 
features of Cole-Cole plot resemble the ones observed experimentally in complex systems. Since 
algebraic waiting time distribution of jumps naturally gives rise to Davidson-Cole kind of fre- 
quency dependence of dielectric function, the present formalism should find use in interpretation 
of existing experimental results. Work in this direction is under progress. 
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Fig. 1: a) Schematic illustration of Kubo model. Jumps are restricted on a circle. 

b) Schematic illustration of Ivanov- Anderson model. Jumps are isotropically distributed on a 

sphere. 

Fig. [2t Correlation functions of Ivanov- Anderson model with A = 60° against dimensionless 
time, j t t. Red lines represent C^t) and black lines represent C\(t). a = 0.5 and 0.8 from 
top to bottom. Thin lines indicate the results of the normal diffusion. Dashed lines represent 
asymptotic algebraic time dependence of Eq. (T2oD . 

Fig. [3} Correlation functions of Ivanov- Anderson model with A = 170° against dimensionless 
time, j T t. Red lines represent C^t) and black lines represent C\(t). a = 0.5 and 0.8 from 
top to bottom. Thin lines indicate the results of the normal diffusion. Dashed lines represent 
asymptotic algebraic time dependence of Eq. ( 1261) . Dots represent the results of Mittag-Leffler 
function of Eq. (127) 

Fig. [4J Log linear plot of Fig. 2. 

Fig. [5t Correlation functions of Ivanov- Anderson model for a = 0.8 with various values of 
jump amplitude against dimensionless time, 7 r i. Red lines represent C^if) and black lines and 
symbols represent C\{t). 

Fig. ED Cole-Cole diagram of Ivanov- Anderson model for A = 30°. 

Fig. O Cole-Cole diagram of Ivanov-Anderson model for a = 0.8. Solid lines indicate 
A = 10°, 60°, 90°, 150°, and 180° (almost overlaps with the line of A = 150°) from bottom to 
top at e' = 0.2. Dashed line represents the Debye relaxation of a = 1.0. 
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Fig. [5J Cole-Cole diagram of Ivanov-Anderson model for a = 0.5. Solid lines indicate 
A = 10°, 60°, 90°, 150°, and 180° (almost overlaps with the line of A = 150°) from bottom to 
top. Dashed line represents the Debye relaxation of a = 1.0. 

Fig. [9t Correlation functions of Kubo model with A = 60° against dimensionless time, 7 r i. 
Red lines represent Cz(t) and black lines represent C\(t). a = 0.5,0.8 from top to bottom. 
Thin lines indicate the results of normal diffusion. 

Fig. HOt Correlation functions of Kubo model with A = 170° against dimensionless time, 
7 r t. Red lines represent C^t) and black lines represent Ci(t). a = 0.5, 0.8 from top to bottom. 
Thin lines indicate the results of normal diffusion. 

Fig. lilt Correlation functions of Ivanov-Anderson model for a = 0.8 against dimensionless 
time, t/r a . Large amplitude jump of A a = 60° and small amplitude jump of A& = 5° coexist. 
The frequency of small amplitude jump is 5 times larger than that of large amplitude jump, 
T a/ T b — 5. Red lines represent (^(t) and black lines represent C\{t). Thick lines indicate the 
exact numerical results. Thin lines indicate the approximate results of eqs. f l50|) - fl5T|) . Short- 
dashed lines indicate the result of small amplitude jump alone. Long-dashed lines indicate 
the result of large amplitude jump alone. The dashed-dotted lines are the results assuming 
independent superposition of two processes as explained in the text. 

Fig. 1121 Cole-Cole diagram of Ivanov-Anderson model for a = 0.8. The solid line represents 
the result when large amplitude jump of A a = 60° and small amplitude jump of A5 = 5° 
coexist. The frequency of small amplitude jump is 5 times larger than that of large amplitude 
jump, T a /rb = 5. The dashed line indicates the result when only the large amplitude jump is 
present. 
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Fig. 113b 7~i / r 2 against the amplitude of the large amplitude jumps for the fixed jump amplitude 
A b = 5° of small amplitude jumps. The frequency of small amplitude jump is 5 times larger 
than that of large amplitude jump, T a /r b = 5. The solid line indicates the results of a = 0.8. 
Dots indicate the results of a = 0.5. The red dashed line indicates the results for the normal 
diffusion of a = 1.0. 

Fig. 114b Correlation functions of Ivanov-Anderson model for a = 0.8 against dimensionless 
time, i/r a . Large amplitude jump of A a = 60° and small amplitude jump of A b = 5° coexist. 
The frequency of small amplitude jump is changed. Red lines represent C*2(t) and black lines 
represent Ci(t). T a /r b = 0.1,r a /r b = 0.5, and r a /r b = 1.0 from top to bottom. 
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Fig. 1, K. Seki, B. Bagchi, and M. Tachiya 
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Fig. 2, K. Seki, B. Bagchi, and M. Tachiya 
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Fig. 3, K. Seki, B. Bagchi, and M. Tachiya 
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Fig. 5, K. Seki, B. Bagchi, and M. Tachiya 
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Fig. 6, K. Seki, B. Bagchi, and M. Tachiya 




Fig. 7, K. Seki, B. Bagchi, and M. Tachiya 
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Fig. 9, K. Seki, B. Bagchi, and M. Tachiya 
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Fig. 10, K. Seki, B. Bagchi, and M. Tachiya 




Fig. 11, K. Seki, B. Bagchi, and M. Tachiya 
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Fig. 12, K. Seki, B. Bagchi, and M. Tachiya 




Fig. 13, K. Seki, B. Bagchi, and M. Tachiya 




Fig. 14, K. Seki, B. Bagchi, and M. Tachiya 
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